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^ : Abstract 

G ' Using a simple model, a new sphaleron solution which incorporates finite 

^-, ■ 

fermionic density effects is obtained. The main result is that the height of the 

potential barrier (sphaleron energy) decreases as the fermion density increases. 

^ \ This suggests that the rate of sphaleron-induced transitions increases when 

(<— ^ I the fermionic density increases. However the rate increase is not expected 

T-H ■ 

^D \ to change significantly the predictions from the standard sphaleron-induced 
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baryogenesis scenarios. 
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Nonconservation of baryon number in the standard model through quantum effects 
(anomahes) is well known [Q. While instanton mediated baryon decays are negligible, the 
same cannot be said for transitions occuring because of monopole catalysis 0, high tem- 
perature 03, high densities |||,§], or in the presence of a heavy particle ^ (for an excellent 
review of these four mechanisms, see ref. [^). 

The basis for all these transitions is the level crossing phenomenon [§ which is usually 
illustrated by looking at adiabatic changes in the gauge field configuration and at the accom- 
panying variation in the energy levels for the fermions resulting in a change in the fermion 
number. This description neglects the effect of the fermion back-reaction: a change in the 
fermion density can introduce a change in the gauge field configuration, just as a change 
of gauge field configuration can change the fermion density. This is most easily seen in the 
Schwinger model where this back-reaction of the fermions is responsible for oscillation in 



the fermion number ||10|. The fermion back- reaction is a purely quantum mechanical effect 
being a direct consequence of the anomaly equation. Since the focus of fermion number vio- 
lation has been in the study of solutions to the classical equation of motion (e.g. instantons 
and sphalerons), it is, therefore, not surprising that little attention has been paid to this 
back-reaction. Further, it might be very difficult to properly take into account the fermion 
back-reaction in realistic 3+1 dimensional theories since the resolution of even seemingly 
straightforward related issues, like the gauge invariance of the free energy at finite temper- 
ature and fermionic density [|lT|, require a careful treatment of non-perturbative effects to 
be properly resolved |]T2|. 



Fortunately, the situation is simpler in 1+1 dimensional models where one can, through 
bosonization |jl3|, take into account the fermion back-reaction at the classical level. The 
present work illustrates some of the non-trivial effects of this fermion back-reaction using 
a simple model which has been extensively studied in the past, namely the Abelian Higgs 
model axially coupled to fermions (see for example PJI^-p!8|) . 

The Lagrangian density describing this model is 



£ = i,ir{d, - %ei''A,)i, - -F,,F^^ + (D^0)(D»* - A(|0|2 - c'f (1) 

where D^ = 9^ — ieA^ and the space coordinate extends from — L to L. This model is to be 
regulated such that the gauged current, ip'j'^'j^ip is conserved while the vector current obeys 
the anomaly equation 

d^i,Y4^^d,J>^ = -^e'^''F,,. (2) 

Instead of working directly with the fermionic Lagrangian, it is preferable to use the 
Bose-equivalent form ||T0|Jl3|Jl6| : 



^ = 1 (d,X - j^A,) - \f,^F^'^ + (D,0)(Z^»* - A(|0p - c'f (3) 



where the mass term for the photon has to be included in order to give the bosonized 
Lagrangian the correct symmetry. The properly regularized vector current is then 



which obeys the anomaly equation (^. 

The main reason for using the Bose-equivalent formulation is that the anomaly equation 
of the fermionic theory is present at the classical level, provided the Lagrangian is properly 
regularized. This ensures that the level-crossing effects connected to the anomaly are in- 
cluded in the classical equations of motion. This is preferable to the usual analyses (e.g. see 



T7|-|19[) where fermions are essentially ignored except for the fact that, through the anomaly 
equation, a change in the Chern-Simons number of a given gauge field configuration is known 
to be accompanied by a corresponding change in the fermion number. 



The static energy density is easily obtained in the x = ^o = gauge [|T6 
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KW-c'Y. (5) 



As a result of this gauge choice, the fermion number and the Chern-Simons number are 
identical. 



Suppose one ignores the explicit mass term for the photon, which is equivalent to ignoring 
the contribution from the fermions. One then finds that £ has an infinite number of local 
minima such that 

-kNcs 



^1 



eL 

cexp(iey4ix) (6) 



where Ncsi an integer, is the Chern-Simons number. Following [[I7i|l8| , one can construct a 



noncontractible path interpolating between two vacuum states pSj using the parametriza- 
tion: 

(P{x;t) = cexp{ieAix) [cositt + isimfr^^z)] (7) 

where z = a/Ac^ sin(7rr) x. As explained by Manton ||2^ and Carson |T^, the sphaleron 



configuration is obtained through a minimax procedure as follows: the set of path 
{(j){x; t) , Ai{t)} are finite-energy field configurations that interpolate between two vacuum 
states as r runs from to 1. Therefore, as a function of r, there must exist a point along the 
path where the energy reaches a maximum, -Emax- By considering the set of all such paths, 
one can find a function ^{z) for which the energy -Emax is minimal. This configuration is the 
desired sphaleron which is a solution to the static equations of motion that corresponds to 
a saddle point of the energy functional for this system. The importance of these sphaleron 
configurations is that they are the main contributor to baryon-number violating processes 
occuring at finite temperature 0. 

It is straightforward to show that the sphaleron configuration for the Abelian Higgs 
model without fermions is given by ^{z) = tanh(2;) and r = ^ in the limit as L — >■ cxd. The 
energy along the corresponding path is then 



E(t) 



sVx, 
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sin'^ TIT 



(8) 
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This results in a periodic effective potential as function of the fermionic density (Chern 
Simons number) as shown by the dashed line in figure 1. 



In order to consider the effects of including the fermions in the system, one may reintro- 
duce the mass term for the photon. Naively, the corresponding effective potential is then as 
illustrated by the solid line in figure 1 p0|-[2^, where the periodic potential obtained before 
is simply added to a pure fermionic contribution. This results in potential barriers of varying 
height separating the various distinct local minima up to a maximum value for the fermion 
density. The computation of transition rates between a state having N fermions to a state 
having A^ + 1 fermions is now more complicated since the transition is between two states 
having different energies rather than a true vacuum-to-vacuum transition. 

Rather than computing the transition rate as a function of the fermionic density, suppose 
one simply wants to compute the critical density. Since the relevant states have a non- 
vanishing energy, it is useful, at finite density and in the L ^ oo limit, to subtract an 
infinite constant from the energy and consider instead 

/oo 
dx{£{N,T)-S{N,T = 0)) (9) 

-oo 

Because the new effective potential was obtained by assuming that one could simply add 
two contributions, it may seem reasonable to assume that the same function ^{z) = tanh(2;) 
found before still corresponds to the sphaleron configuration. Thus, one finds 



8v^, 



,3 



E{t) = 27rrn H sin"* vrr (10) 

o 

where n = N/L is the finite fermion density taken to be finite even in the L ^ oo limit. 
This result shows that it is possible to find a local extremum of E{t) only if Ai does not 
exceed a critical value given by 327rAc^/3e-\/3. 

However, this result indicates that such a naive way of deriving an effective potential at 
finite density is not correct since, as has been previously found P, p3| , p!6| , no local minima 



of the energy exists when Ai exceeds 2\^7i\c^ /3e\^. Furthermore, this derivation doesn't 
take into account the fermion back- reaction and the function ^{z) is not a static solution 
to the complete set of equations of motion of the original Lagrangian (eq. ^). Writing 
= $exp(ip), the static equations of motion are: 



eAi - 27r<l>^(9^p - eAi) = (11) 

d^[^\d^p-eAi)] =0 (12) 

92$ - $(9^p - eAif - 2A$($2 - c^) = (13) 

Any solution to these equations has to be such that Ai is spatially constant and that 

The periodicity requirement on p then implies that the various local minima of the energy 
functional are such that 

Ncs = ^AL = ^""^ (15) 



L \ 27r<l>2 



is no longer an integer. In other words, the minima no longer coincide with pure gauge 
configurations. Furthermore, the difference between two adjacent minima of this quantity 
is also different from unity. This complicates the search for a non-contractible path joining 
two adjacent minima since it becomes very difficult to find a parametrization consistent 
with the periodic boundary condition. However, a static solution to the equations of motion 
resembling the Abelian Higgs sphaleron solution can still be found. If one combines equations 
(0) and (p!3D , one gets 



dx^ 7r2<|)3 



2A$($^ - c') (16) 



which, in the L ^ 00 limit, has a single non- homogeneous solution [|J] as well as two 
homogeneous solutions. To see this, it is convenient to parametrize Ai in terms of a new 
variable 7 such that 



--(-i)#^ • ^<-^<-^ (-' 

This parametrization is one-to-one on the defined interval and such that y4i(7 = 0) =0 
while ^1(7 = 1) is the critical density above which no local minima of the energy exists. 
With this parametrization, the inhomogeneous solution (sphaleron) is given by 
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^sph = cJy + (1 - 7) tanh^ ( ^A(l - 7)cx ) (18) 



and the two homogeneous solutions are 



$i = cjl-| 



$2 = —^7 + 73^^4^ (19) 



c 
76' 



The difference in energy between the sphaleron solution and the lowest energy homogeneous 
solution ($1) is 8c^wA(l — 7)/3 which correctly reproduces the Abelian Higgs result when 



the fermion density is zero (7 = 0), vanishes at the critical density (7 = 1), and remains finite 
for all other values of 7. The result for the effective potential is given in Figure 2. There are 
three important characteristics of this effective potential; firstly, its overall positive curvature 
due to the finite fermionic density; secondly the existence of a critical density above which 
no stable solution can be found; and lastly, and perhaps most importantly, the decrease of 
the sphaleron energy (the barrier height) with increasing fermion number. Of course, the 
important question is the extrapolation of these results for 3+1 dimensional theories. Firstly, 
in order to compute transition rates at finite temperature, the relevant quantity is the free 
energy which also has a generally positive curvature as a function of the fermion density 
Hj20|-p2|,p5i . Secondly, in 3+1 dimensional models there also exists a critical density [P,|l, 



just as in the 1 + 1 dimensional model discussed above. Hence, it is probably safe to assume 
that the height of the barrier between local minima (the energy of the sphaleron) decreases 
as the fermionic density increases so that it vanishes at and above the critical density. 
From this it follows that, as the fermionic density increases, sphaleron induced transitions 
would occur even more rapidly than the rate suggested by conventional calculations, with 
the finite density effects becoming the dominant factor when the fermionic density becomes 
comparable to or greater than the critical density. The critical density computed in the 
standard model |^,^ is about 12 orders of magnitude greater than the nuclear density which 
is comparable to the total energy density close to the electroweak phase transition in the 
standard Big Bang model [^. However, the critical density is the net fermionic density 



(i.e. ferinion minus anti-fermion) which is hkely to be much smaller than the overall density 
given the small baryon to photon ratio obtained from primordial nucleosynthesis. Thus it 
appears that the rate increase of sphaleron transitions due to finite density effects would be 
negligible. Nonetheless, it might still be desirable to compute this rate increase precisely 
and to do so it will be necessary to know in what way the standard sphaleron ^^ is modified 
at finite density. It is likely that non-perturbative effects would play a crucial role thereby 
leading to an impossibly complicated analytical solution. However, a numerical solution 
should be attainable. 
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FIGURES 
FIG. 1. Schematic representation of the naive effective potential as a function of the fermion 

density. The dashed line represents the abelian Higgs contribution which has an infinite number of 

local minima. The dotted line is the pure fermion contribution. The solid line, obtained by adding 

the abelian Higgs sphaleron contribution to the pure fermion contribution is not the true effective 

potential. 

FIG. 2. Schematic representation of the real effective potential as a function of the fermion 
density. The main thing to note is that the height of the potential barrier decreases with increasing 
fermion density until it vanishes at the critical density. 
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This figure "figl-l.png" is available in "png" format from: 
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